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C/3 ! Abstract 
£h 

The theory of linear wave structures generated in Bose-Einstein condensate 
flow past an obstacle is developed. The shape of wave crests and dependence 
of amplitude on coordinates far enough from the obstacle are calculated. The 
results are in good agreement with the results of numerical simulations obtained 
earlier. The theory gives a qualitative description of experiments with Bose- 
Einstein condensate flow past an obstacle after condensate's release from a trap. 
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: 1 Introduction 

in 

After Bose-Einstein condensation of dilute gases was obtained experimentally, great 
O ! attention has been devoted to excitations existing in such systems (see, e.g., pQ). For 

example, it is well known that sound waves can propagate through the condensate 
with repulsive interaction between atoms. Such a condensate can occupy considerable 
volume inside a trap, therefore, if a wave length is much smaller than the condensate 
size, one can consider condensate ad locally homogeneous one with a uniform density 
n Q in an undisturbed state. As was shown by Bogoliubov [2] , the linear waves have the 
dispersion law 



2m 



where 
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is a speed of sound in a long wave limit k — > 0, and g denotes the effective coupling 
constant 

47ih 2 a s 

9 = 3 

m 

arising due to s-wave scattering (with the scattering length a s ) of atoms with mass m. 
The repulsion of atoms corresponds to g > and a real speed of sound c s . For small 
wave numbers the dispersion law (CQ) corresponds to the sound waves propagating with 
speed c s , 

to ^ c s k, k-^0, (4) 
but for large k it gives the quantum dispersion law of free particles with mass m, 

hk 2 



2m 



oo. (5) 



Transition from one limiting case to another takes place at k ~ £ 1 where £ is a 
characteristic parameter called a healing length 

{=-»=* (6) 
y/2mgn V2mc s 

If the amplitude is not small and the condensate is not homogeneous, in the mean- 
field approach dynamics of the condensate is described by the Gross-Pitaevskii equation 

where is a "condensate wave function", and V(r) is an external potential (for 
example, it might be the confining potential of a trap). In the region where condensate 
is homogeneous, potential V is constant; hence, it can be excluded from the equation 
by multiplication of ^-function by an inessential phase factor. Linearizing equation 
(J7J) with respect to a state with a constant density n = \ip\ 2 yields solutions in the 
form of plane waves with the dispersion law ([Tj). But appearance of the nonlinear 
term in (jTj) leads to new kinds of nonlinear excitations — vortices and dark solitons. 
Linear sound waves of the described origin were observed in the experiment [3] and 
dark solitons in jU [5]. If we consider evolution of the condensate with effectively one 
spatial dimension and the initial disturbance of density is large enough, then the pulse 
of rarefication will decay into a series of dark solitons [5]. But if the initial disturbance 
of density is positive, it will lead to wave breaking of the pulse with formation of the 
dispersive shock wave [3 [8] , as it has been found experimentally in [9J [10] . 

Linear waves, dark solitons and shock waves referred above can be described in 
the first approximation as waves depending only on time and one space coordinate. 
But situation becomes completely different if the wave behavior depends on two space 
coordinates. In particular, the experiment [TT] (see also [12]) has been made where 
condensate released from a trap interacted during its expansion with a small obstacle 
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generated by a laser beam. As was shown experimentally, the flow past an obstacle 
generates wave structures in the form of stationary "moustaches" attached to the ob- 
stacle. Theoretical analysis [321 [H] shows that the wave structure generated in the 
flow can be subdivided into two regions separated by the Mach cone corresponding to 
the sound velocity c s in the long wave limit, i.e. by lines drawn at the angle 9 to the 
direction of the flow so that 

u M 

where M is the Mach number. As was shown in [T3J [H] , inside the Mach cone dark 
solitons are generated. Their slope depends on their depth — the most shallow solitons 
lie near the Mach cone and the deepest ones are located near the flow axis going 
through the obstacle. In the paper [13] the exact stationary solution of the Gross- 
Pitaevskii equation was found and it was shown by numerical calculation that this 
solution describes dark solitons generated inside the Mach cone. Linear waves are 
characterized by the dispersion law (CQ), where k means the module of the wave number. 
This dispersion law shows that both phase and group velocities are greater than the 
sound velocity c s in a long wavelength limit. It means that linear wave structure formed 
by a wave packet can only occupy the region outside the Mach cone defined by the 
equation ([8]). For the first time these waves described by the Gross-Pitaevskii equation 
were observed in the numerical experiment [To] . In the paper [12] analogous numerical 
simulation was carried out for description of the structures observed in the experiments 
with the condensate flow past obstacle. Some of the elementary characteristics of these 
structures were derived from the dispersion law ([1]). In the paper [16] it was noticed that 
these waves have a physical origin analogous to Kelvin's ship waves [17] generated by 
a ship moving in a still water. It means that similar computation methods [TB], [TH] can 
be used in the case of Bose-Einstein condensate. In [TH] such a theory was developed 
and shape of waves crests was obtained. It was shown that this theory agrees very well 
with numerical simulations. But it cannot give the dependence of the wave amplitudes 
on space coordinates. The aim of this paper is to develop a complete theory which 
can give not only geometrical but also dynamical properties of the wave structures 
generated outside the Mach cone by the condensate flow past obstacle. 



2 The theory of stationary wave structures gener- 
ated by condensate flow past obstacle 

As a basis of the theory, we take Kelvin's method in the form presented in [18J. Ac- 
cording to this theory, the wave structures result from interference of linear waves 
emitted by the source at all previous moments of time. We take a source (obstacle) 
at rest, in agreement with the experiment [TT], and suppose that condensate moves 
with a constant speed Uq in the direction of the x axis. Since a stationary flow cor- 
responds to the potential V(r) turned on in the infinite past, one can consider this 
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physical condition by the assumption that the potential rises in time according to the 
law V(r,t) = exp(et)V(r). Afterwards one should take the limit e — > 0. 

It is convenient to transform the Gross-Pitaevskii equation if)(r,t) to a "hydrody- 
namic" form by means of the substitution 

V>(r,t) = y/n(T,t)exp U j U(r',t)dr') , (9) 

where n(r, t) is a local density of atoms in condensate, U(r, t) is a velocity field of 
condensate flow. To simplify the notation, let us introduce non-dimensional variables 

r ; t, h = —, U = -, (10) 



V2^ V2£ ' n 

and dimensionless potential of the obstacle 

V-(r) = -Lv{V2£r). (11) 

where no is a density of undisturbed condensate (infinitely far from the obstacle). 
Upon substitution of (jHD, ( JTOl and (1TTI) into ((Zj) and separation of real and imaginary 
parts, we obtain equations for the condensate density n and velocity components U = 
(u(x,y,t),v(x,y,t)) in a hydrodynamic form: 

n t + (nu) x + (nv) y = 0, 

u t + uu x + vu y + n x + {^±pL _ ^±%) ^ = -y„ e - (i2) 

/ n 2 x + n^j, + n vv \ T r 
v t + + w„ + n y + ^ g ^ 2 • ^-^ J = -V„e , 

where tildes are omitted for convenience of the notation and the indices stand for 
derivatives with respect to appropriate variables. 

We are interested in linear waves propagating in a uniform flow with n = 1, 
u = M, Vq — 0. Therefore we introduce new variables 

n = 1 + ni, u — M + Ui, v — V\ (13) 

and linearize the system (ITS]) with respect to small deviations ni, Ui, V\ from a uniform 
flow. As a result we arrive at a linear system 

ni,t + u 1)X + Mni tX + v hy = 0, 
u 1;t + Mu 1;X + n liX - l{n 1>xxx + n liXyy ) = -V x e £t , (14) 
v x ,t + Mvi jX + n l)V - l(nx )Xxy + n Xmy ) = -V y e £t . 
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Figure 1: Polar coordinates defining the radius vector r and the wave vector k normal 
to the wave front which is shown schematically by a curved line. 



A stationary solution corresponds to the limit e — > 0, where all the variables depend 
on time as exp(ei). Fourier transform for spatial variables yields the solution 

, , k 2 <$>(k x ,k y )e< k * x+k ^ dk x dk y 

u± — —c 

where 

-jkr 



' ' (e + ik x M) 2 + k 2 (l + k 2 /4) (2vr) 2 ' 1 ; 



$(k) = / V(r)e- tkr dr. (16) 

Far enough from the obstacle its potential can be considered as a point-like one, 
V(r) = VqS(t); hence we have $(k) = Vq. It is convenient to introduce polar coordi- 
nates (see Fig. 1) for vectors r and k: 

x = rcosx, 2/ = r sinx; k x = —kcosr], k y = ksmr]. (17) 

Then simple transformation of equation (fl5|) gives 

V £t r f°° ke- ikrcoB ^+^dkdrj , 1C . 



m = — e 

IT 



St 

-7T Jo 



k 2 — k 2 — i5 cos r) 
where 5 = 8Me/k is a small positive quantity and 



k = 2^/M 2 cos 2 ?/ - 1. (19) 

One can divide the integration interval over r\ into two parts: f dr) — J^ 2 dr]+ drj. 
After replacement r( = 77 — 7r in the second integral, one can notice that the integrand 
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function converts to its complex conjugate one. Hence we can write the integral (fT8l) 
in the form 

..^frf^, (20) 



7T 
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k 2 — k\ — i5 cos 77 



Integration over should be taken along the positive axis. The pole 

. 5 cos 77 



fc = V fc o + ^ cos 77 = k +i - , (21) 
v 2fc 

lies in vicinity of a real axis, so we can expect that it gives the main contribution 
into the integral. This can be proved by the following argumentation. We add to 
the integration contour along real positive k axis also integration along the positive 
or negative imaginary axis and a quarter of a circle with infinite radius to make the 
contour closed. If cos(x + if) < 0, we use contour along the boundary of the first 
quadrant of complex fc-plane, on the other hand, if cos(x + if) > 0, we use the contour 
along the boundary of the fourth quadrant. In both cases integration along a quarter 
of a circle vanishes when radius \k\ of the circle tends to infinity. Thus, we have 
no contribution of this part of the integration contour. The pole (|2ip gets inside the 
integration contour only if one take the integral along the boundary of the first quadrant 
what corresponds to the case cos(x + i]) < 0. In this case the result of integration along 
the real positive axis k is equal to contribution of the residue and integration along the 
positive imaginary axis. If we have cos(x + 77) > 0, this integral is equal to the integral 
along the negative imaginary axis only. In both cases an estimation of integral along 
the imaginary axis gives 

ke -krco B ( X +r,) dk l 

oc — , (22) 



/ k 2 + kl r 2 : 

which means that this contribution decreases as r~ 2 with distance r from the obstacle. 
As we shall see below, the contribution of the pole decreases as r -1 / 2 . It means that 
far enough from the obstacle we can neglect the contribution of integration along 
the imaginary axis. Thus, a noticeable wave structure appears in the region where 
cos(x + v) < 0; i- e - the angle \i = n — (x + if) between vectors r and k should be acute. 
The contribution of the pole (1211) is equal to 

2V T f n/2 

m = 



71 



Im f e- ifcrcos(x+r?) d77, (23) 

J-n/2 



where k is determined by the equation ffl9l) (index "0" is dropped here). 

Far from the obstacle where the phase kr = rs is large enough, we have large values 

of 

s(r]) = k(r])cos(x + v), ( 24 ) 

and integral (1231) can be estimated by a standard method of stationary phase. Condition 
ds/dr] = gives the equation for the point of stationary phase and it can be easily 



6 



transformed to the form 

2M 



tan// = - sin2r/ (25) 



k 

or, with the help of equation for /x, it gives the expression for \ 



(1 + Ar72)tan77 
tmX= jU-(l + *»/2) - (26) 

Taking into account equation f[T9"l) , we find 

k . , 

COS ^ = 2[(M 2 -2)fc 2 + 4(M 2 -l)]V2- (27) 

With account of (1251) . we get the expression for a second derivative of the phase 

^± = 8^[(M 2 - 2)A; 2 + 6(M 2 - 1)]. (28) 

As a result, the expression for the condensate density f[2"3"j) takes the form 

/2l : [(M 2 -2)A; 2 + 4(M 2 -l)] 1 / 4 / tt\ 

7ll = VoA/ FTTTn 777; TTT^ -77777 COS I KT COS /i . (29) 

V txt [(M 2 -2)P + 6(M 2 - I)} 1 / 2 V p 4/ y ' 

First of all, let us check that the results obtained above agree with the theory 
developed in [16J. Since the wave crest corresponds to the constant value of the phase 
4> = kr cos /i, we can find the equation for the absolute value of the radius- vector 

[(M 2 - 2)A; 2 + 4(M 2 - l)] 1 / 2 . (30) 



k cos fi k 3 

Calculating cosx and sin% with the help of equation ( 1261) . we obtain equations for lines 
of stationary phase coinciding with those found by another method in |16j : 

40 

x = r cosx = 77 cos 77(1 — M 2 cos 2rj), 

u < 31 > 

y = r sin v = — sin r](2M 2 cos 2 77 — 1). 
k 

Since the wave vector k is determined as a function of 77 by equation (Tl9|) . equations (j3~Ij) 
give the wave crests shape in a parametric form where parameter 77 varies in the interval 

- arccos(l/M) < 77 < arccos(l/M). (32) 

For small 77 we have 

* (2M 2 - 1)0 ^ 2 



2(M 2 - I) 3 / 2 
7 



+ -TTTTo T^V 



2VM 2 - 1 4(M 2 -1) 3 / 2 
(2M 2 - 1 



2 , , ( 33 ) 
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Figure 2: Wave structure calculated according to equation (1291) with M = 4; it is 
generated outside the Mach cone and only its upper part with y > is shown. 



that is the lines of stationary phase take here a parabolic form 

^-27fe + W^V- (34) 

The limiting values rj = iarccos (1/M) correspond to the lines 

- = ±VM 2 - 1, (35) 

y 

i.e. far from the obstacle wave crest curves approach to the straight lines parallel to 
those forming the Mach cone ([8]). 

An example of the wave structure determined by equation ff29l) with M = 4 is 
shown on Fig. 2. 

If M > 1, there exists a region where M cos 77 3> 1; all the equations get here a 
simple form. In particular, we have x — 71 ~ 2r7, /i = rj = (tt — x)/2, so that 

A; = 2Msin|>l (36) 
and as a result the expression for condensate density (T29]) takes the form 

^^ y °\^ COS ( 2Mrsin2 f-l)- ^ 
8 
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Figure 3: Oscillations of the condensate density in a stationary wave in front of the 
obstacle: x < 0, y — 0. The condensate flow corresponds to the Mach number M = 4. 



The expressions f l3T|) take here a parabolic form 

M 5 

x = ~7777 + 777?/ 



(3? 



2M 20' 

Naturally, this result can be obtained from (l32l) in the limiting case M ^> 1. 

In the region in front of the obstacle where y — 0, x < 0, the perturbations of the 
condensate density take the simplest form. Here we have 



k = 2VM 2 - 1, 
i.e. the wave length A = 2n/k is constant and 



( M 2 _ !)l/2 
T2i = 2Vo4 / — — COS 



-2VM 2 -lx - 



y = 0, x < 0. 



(39) 



(40) 



The plot illustrating this dependence is shown in Fig. 3. All the results obtained above 
are in a good agreement with a numerical simulations carried out in [16] . 



3 Conclusions 

The theory developed here is valid as long as the amplitude of the wave structure is 
much smaller than unity (in dimensionless variables), i.e. 

r > V 2 /M. (41) 
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Far enough from the obstacle this condition is always fulfilled. But for comparison of 
our theory with experiment it is also necessary that the condensate flow be uniform 
in the region under consideration. The estimates made in [T3j show that one can 
satisfy this condition for a free flow of condensate released from the trap by locating 
the obstacle far enough from the center of the trap and sufficiently long time after 
condensate's release. As far as we can see, this conditions were not fulfilled in the 
experiment [HI [T2]. In particular, the obstacle was situated too close to the center 
of the flow and its size was too large, so that the shadow region appeared behind 
it which was not filled by the condensate. We believe that for this reason no dark 
solitons predicted in [13J were observed experimentally. Nevertheless the wave structure 
in [11] [12] looks similar to one shown in Fig. 2. Thus we believe that the theory 
developed above gives at least qualitative description of the wave structures observed 
experimentally. 
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